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We report in this paper the proofs that the pulse shape analysis can be used in some bolometers to identify the nature of the 
interacting particle. Indeed, while detailed analyses of the signal time development in purely thermal detectors have not produced 
so far interesting results, similar analyses on bolometers built with scintillating crystals seem to show that it is possible to distinguish 
between an electron or y-ray and an a particle interaction. This information can be used to eliminate background events from the 
recorded data in many rare process studies, especially Neutrinoless Double Beta decay search. Results of pulse shape analysis 
of signals from a number of bolometers with absorbers of different composition (CaMo04, ZnMo04, MgMo04 and ZnSe) are 
presented and the pulse shape discrimination capability of such detectors is discussed. 
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1. Rare event searches 

Rare event studies, such as the search for Neutrinoless Dou- 
ble Beta decay (J3/3(0v)) [ 1 1 or the identification of Weakly In- 
teracting Massive Particle (WIMP) interactions with ordinary 
matter ||2), are of extreme interest in astroparticle physics, since 
they would imply new physics beyond the Standard Model. In 
both cases, as in all the rare event studies, spurious events are 
a limiting factor to the reachable sensitivity of the experiment. 
Unfortunately natural radioactive background is often present 
in the detector itself or in the materials surrounding it, no mat- 
ter how much one can try to reduce it with shieldings, selection 
of materials and complicated purification techniques. In order 
to handle the residual unavoidable background, all the envis- 
aged approaches require both a good energy resolution (which 
always helps in the comprehension of the different structures 
of an energy spectrum) and the capability to identify the nature 
of the projectile that interacted with the detector. Indeed, the 
searched event has always well defined signatures helping to 
distinguish it from background, for instance two electrons with 
a fixed sum energy in the case of the /3/3(0v) . 

Bolometers [3| are based on the detection of phonons pro- 
duced after an energy release by an interacting particle and can 
have both an excellent energy resolution and extremely low en- 
ergy threshold with respect to conventional detectors. They can 
be fabricated from a wide variety of materials, provided they 
have a low enough heat capacity at low temperatures, which 
is the only requirement really unavoidable to build a working 
bolometer. The latter is a priceless feature for experiments that 
aim at detectors containing particular atomic or nuclear species 



to optimize the detection efficiency. If other excitations (such 
as ionization charge carriers or scintillation photons) are col- 
lected in addition to phonons, bolometers have already shown 
to be able to discriminate nuclear recoils from electron recoils, 
or a particles from /3 particles and y-rays. In this paper we will 
report on the possibility to obtain similar results just by pulse 
shape analysis, without the requirement of a double readout for 
phonons and ionization or scintillation light. 

2. Bolometric Technique and Scintillating Bolometers 

Bolometers can be essentially sketched as a two-component 
object: an energy absorber in which the energy deposited by a 
particle is converted into phonons, and a sensor that converts 
thermal excitations into a readable signal. The absorber must 
be coupled to a constant temperature bath by means of a weak 
thermal conductance. 

Denoting by C the heat capacity of the bolometer, the temper- 
ature variation induced by an energy release E in the absorber 
can be written as 



AT = 



C 



(1) 



The accumulated heat flows then to the heat sink through the 
thermal link and the absorber returns to the base temperature 
with a time constant r = C/G, where G is the thermal conduc- 
tance of the link: 



E 

AT(t) = -e 



(2) 
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In order to obtain a measurable temperature rise the heat ca- 
pacity of the absorber must be very small: this is the reason 
why bolometers need to be operated at cryogenic temperatures 
(of the order of 10-100 mK). 
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A real bolometer is somewhat more complicated than the 
naive description presented above. It is made of different el- 
ements and it is therefore represented by more than one heat 
capacity and heat conductance. As such, the time development 
of the thermal pulse is characterized by various time constants. 
In principle, if the bolometer performs as an ideal calorime- 
ter and if the conversion of the energy into heat deposited by 
the particle is instantaneous (as assumed in equation [2]), then 
the device is insensitive to the nature of the interacting parti- 
cle. Although this situation is generally very far from reality, it 
is however true that the small differences are difficult to detect 
and the goal has been so far achieved only relying on compli- 
cated solutions. Among these are scintillating bolometers. 

The concept of a scintillating bolometer is very simple: a 
bolometer coupled to a light detector H. The first must consist 
of a scintillating absorber thermally linked to a phonon sensor 
while the latter can be any device able to measure the emitted 
photons. The driving idea of this hybrid detector is to combine 
the two available informations, the heat and the scintillation 
light, to distinguish the nature of the interacting particles, ex- 
ploiting the different scintillation yield of /3/y, a and neutrons. 
Dark Matter as well as /3/3(Qv) searches can benefit of this capa- 
bility of tagging the different particles, and more generally this 
technique can be exploited in any research where background 
suppression or identification is important. 

Dark matter experiments look for a very rare signal gener- 
ally hidden in a huge background. The signal is a nuclear recoil 
with an energy of few keV (or less) induced by the scattering 
of a WIMP off a target nucleus. Experiments like CDMS 0, 
Edelweiss [6| or CRESST [7| clearly show that in such energy 
region the background is dominated by /3/y interactions. A sec- 
ond source of background are a decays, contributing through 
energy degraded a's and nuclear recoils. The capability to dis- 
tinguish a nuclear recoil - candidate for a WIMP interaction - 
from a or /3/y clearly allows to improve drastically the experi- 
mental sensitivity. 

A similar approach was proposed also for applications in 
/3/3(0v) searches |4|. More recently, such a possibility has been 
demonstrated to be viable for a number of candidate nuclei 1 8 1 . 
In this case the major interest is the identification of a interac- 
tions. Indeed the other important source of background, namely 
y-rays, is virtually indistinguishable from the /3/3(0v) signal. 
The suggested way to eliminate the problem of y-rays con- 
tribution is to study isotopes with a transition energy above 
2615 keV. This corresponds in fact to the highest energy y-ray 
line from natural radioactivity and is due to 208 T1. Above this 
energy there are only extremely rare high energy y's from 214 Bi 
(all the (3/3 active isotopes with Qpp >2615 keV are listed in 
Tab. [T}. Once y-rays are no more a worrisome source of back- 
ground, what is left - on the side of radioactivity- are a emis- 
sions. Indeed a surface contaminations not only can represent 
the dominant background source for /3/3(Qv) searches based on 
high transition energy isotopes, but they are already recognized 
as the most relevant background source in the bolometric ex- 
periment CUORICINO 10 [TUl [TT] and as a limiting factor for 
the experiment CUORE lfT0l[T2l . Both the experiments search 
for the /3/3(0v) of 13l) Te whose transition energy is at 2527 keV, 



Table 1 : Double beta decay isotopes with endpoint energies above the 208 T1 
line. 



Isotope 


Q«s [MeV] 


natural abundance 


116 Cd 


2.80 


7.5 % 


82 Se 


3.00 


9.2 % 


I00 Mo 


3.03 


9.6 % 


96 Zr 


3.35 


2.8 % 


iso Nd 


3.37 


5.6 % 


48 Ca 


4.27 


0.19 % 



therefore in a region where y background (mainly due to Comp- 
ton events produced by 2615 keV photons) can be still impor- 
tant. 

The a contribution to the background in the /3/3(Qv) region 
(i.e. at about 3 MeV) is the following. In the natural chains 
we have various nuclei that decay emitting an a particle with 
an energy between 4 and 8 MeV, their energy is quite higher 
than most /3/3{Qv) Q-values. However, if the radioactive nucleus 
is located at a depth of a few yum inside a material facing the 
detector, the a particle looses a fraction of its energy before 
reaching the detector and its energy spectrum looks as a contin- 
uum between and 4-8 MeV 1131 . A similar mechanism holds 
in the case of surface contaminations on the bolometer. 

This radioactive source plays a role in almost all detectors 
but it turns out to be particularly dangerous for fully active de- 
tectors, as in the case of bolometers. It can be efficiently iden- 
tified and removed with active background suppression tech- 
nique such as that conceived with scintillating bolometers. 

The development of a hybrid detector, able to discriminate 
a particles and optimized for /3/3(Ov) searches, was the main 
purpose of our studies on scintillating bolometers. We tested 
several devices, differing mainly in the scintillating crystal ma- 
terial and size, to study their thermal response, light yield 
and radio purity. The results obtained so far are reported in 
lfT4l[T5l[T6l[T7ll . On our way, we discovered an extremely inter- 
esting feature of some of the tested crystals: the different pulse 
shape of the thermal signals produced by a's and by /3/y's. This 
feature opens the possibility of realizing a bolometric experi- 
ment that can discriminate among different particles, without 
the need of a light detector coupled to each bolometer. In the 
case of a huge, multi-detector array, such as CUORE [12J and 
EURECA 1 18 1, the benefits of employing this technique would 
be impressive: 

• more ease during assembly because the single element of 
the array would be a quite simpler device. 

• fewer readout channels, with not only an evident reduc- 
tion of cost and work, but also a cryogenic benefit (in a 
cryogenic experiment particular care should be devoted to 
reduce any thermal link between room temperature and the 
bolometers working at few mK: the heat load of the read- 
out channels must be taken into account and their reduc- 
tion is always a good solution). 
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• a significant cost reduction, saving money and work that 
would be necessary for the light detectors procurement and 
optimization. 

• no need of light collectors, this would simplify the struc- 
ture of the assembly and it would allow the use of coinci- 
dences between facing crystals to further reduce the back- 
ground. 

As a final remark, it is worth to be mentioned that these de- 
vices could be used also for the measurement of a emissions 
from surfaces, when extremely low counting rates are needed. 
Indeed, due to their lack of a dead layer and their high energy 
resolution, bolometers have an extraordinary sensitivity to low 
range particles like cr's. However, a conventional bolometer 
cannot distinguish the nature of the interacting particle. It pro- 
vides therefore only a limited diagnostic power (especially for 
a particles with energies lower than 2615 keV where the p/y 
induced background dominates the detector counting rate). On 
the other hand, a scintillating bolometer has to be surrounded 
by a reflector to properly collect the scintillation light (therefore 
cannot be faced to a sample whose radioactive emission has to 
be identified). The devices here discussed overcome these two 
difficulties. Traditionally the devices used in this field are Si 
surface barrier detectors. For low counting rates, large area low 
background detectors are needed. Today Si surface barriers de- 
tectors with an active area of about 10 cm 2 , a typical energy 
resolution of about 25-30 keV FWHM, and counting rates of 
the order 0.05 count/h/cm 2 between 3 and 8 MeV are available 
fl9l . A bolometer like those here discussed can easily reach a 
much larger active area, has a typical energy resolution of 10 
keV and a background counting rate in the 3-8 MeV region that 
can be as low as 0.001 count/h/cm 2 . Thanks to the particle iden- 
tification technique discussed in this paper, it can distinguish an 
a emission from a Ply oneQand finally can reject the Ply back- 
ground extending its measurement field to energies by far lower 
than 3 MeV. 

In the following section we report the results obtained with 
the pulse shape analysis on some of the tested crystals. 

3. Detectors, set-up and data analysis 

The results discussed in this paper have been obtained operat- 
ing different scintillating bolometers in an Oxford 200 3 He/ 4 He 
dilution refrigerator located deep underground, in the National 
Laboratory of Gran Sasso (L'Aquila, Italy). The rock overbur- 
den (average depth ~3650 m.w.e. Il20ll ) ensures a strong sup- 
pression of cosmic rays that in our case is mandatory to be able 
to operate the detectors without an overwhelming pile-up. A 
detailed description of the experimental setup can be found in 
ETII . In order to study the pulse shape characteristics of dif- 
ferent materials (in particular those of interest for PP(0v) ), we 
operated a number of scintillating bolometers differing for size, 



This feature is extremely important when the sample which is investigated 
produces a continuous counting rate which could be due either to a emissions 
from a thick contamination or to a ft/y continuum. 



geometry and, of course, the absorbing material. As light detec- 
tor we have used a second bolometer able to absorb scintillation 
photons converting their energy into heat. This was realised us- 
ing as absorber Ge wafers of about 5 g, covered - on the side 
facing the scintillating crystal - with a 600 A thick layer of S1O2 
in order to increase the light absorption. In this way they pro- 
vided measurable thermal signals over an extremely large band 
of scintillation wavelengths. 

Both the scintillating crystal and the Ge wafer were equipped 
with a Neutron Transmutation Doped Ge thermistor (NTD) 
[22], glued on the crystal surface and used as a thermometer to 
measure the heat or light signal produced by particles traversing 
the scintillating crystal. 

A silicon resistance, glued on the crystals, was used to pro- 
duce a calibrated heat pulse in order to monitor the thermal gain 
of the bolometer. This is indeed subject to variation upon tem- 
perature drifts of the cryostat that can spoil the energy resolu- 
tion. In most cases this temperature drift could be re-corrected 
off-line on the basis of the measured thermal gain variation 

£3 EH. 

3.1. Read-out and DAQ 

The read-out ll25ll of the thermistor was performed via a 
preamplifier stage, a second stage of amplification and an 
antialiasing filter (a 6 pole roll-off active Bessel filter 120 
db/decade [26 1) located in a small Faraday cage. The ADC was 
a NI USB-6225 device (16 bit 40 differential input channels). 
For each triggered signal the entire waveform (raw-pulse) is 
sampled, digitized and acquired for the off-line analysis. Since 
all the relevant parameters (including the amplitude) of the trig- 
gered signals are evaluated off-line, a particular care has to be 
dedicated to the optimization of the signal filtering and digiti- 
zation. In the case of the scintillating bolometer the large heat 
capacity of the absorber, coupled to the finite conductance of 
the crystal-glue-thermistor interface results in quite slow sig- 
nals, characterised by a rise-time of the order of few ms^jand a 
decay-time of hundreds of ms (determined by the crystal heat 
capacity and by its thermal conductance toward the heat sink). 
Consequently the sampling rate typically used for the signal is 
1-4 kHz, over a time window of 200-2000 ms. The Bessel fil- 
ter acts mainly as antialising, to avoid spurious contributions 
in the sampled signal. Generally it is preferred to fix its cut- 
off frequency at the lowest value that does not deteriorate the 
signal to noise ratio (i.e. to obtain the best results in terms of 
energy resolution). This results to be a frequency of the order 
of 10 Hz, which is by far lower than what needed for antialias- 
ing purposes. In the studies here presented the Bessel cut-off 
frequency was fixed at 120 Hz in order to exploit the maximum 
available information in the signal bandwidth. 

3.2. Analysis techniques 

Off-line analysis aims at determining the pulse amplitude and 
energy together with several pulse shape parameters associated 



2 The minimum observable signal rise-time is limited by the integration on 
the parasitic capacitance of the signal wires that connect the NTD thermistor to 
the front-end electronics. 
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with each raw-pulse waveform recorded by the data acquisition 
system. Starting from these quantities the physical informations 
that are relevant for the scientific goals can be extracted. 

The first step of the analysis consists in the correct evalua- 
tion of the pulse amplitude. Since thermal pulses are super- 
imposed to stochastic noise, a simple maximum-minimum al- 
gorithm would not give the better achievable resolution. We 
therefore use the Optimum Filter (OF) approach [27]. This al- 
gorithm has proven to provide the best estimate of the pulse 
amplitude and, as a consequence, the best energy resolution. 
The basic concept is to build a filter that, when applied to the 
raw-pulse, produces - as output - a pulse with the best signal to 
noise ratio. The filtered pulse is then used to evaluate the signal 
amplitude. It can be proven that in the frequency domain the 
OF transfer function H( a> ) is given by 



V ' N(oj) 



jujIm 



(3) 



where S(w) is the Fourier transform of the ideal thermal sig- 
nal (reference pulse in the absence of noise), N(w) is the noise 
power spectrum, tu is the delay of the current pulse with respect 
to the reference pulse and K is a proper normalising factor usu- 
ally chosen in order to obtain the correct event energy. 

The role of the optimum filter is to weight the frequency com- 
ponents of the signal in order to suppress those frequencies that 
are more affected by noise. It can be seen from eq. [3] that, in 
order to build the filter, the shape of the reference pulse S(co) 
and the noise power spectrum N(a>) must be known. S(w) is 
usually estimated by averaging a large number of recorded raw- 
pulses, so that the noise associated with each of them averages 
to zero. N(w) is obtained according to the Wiener-Khintchine 
theorem by acquiring many detector baselines in absence of 
thermal pulses and averaging the corresponding noise power 
spectra. 

Once the pulse amplitude has been evaluated, gain instability 
corrections are applied to data. Due to the dependence of the 
detector response on the working temperature, the same amount 
of released energy can produce thermal pulses of different am- 
plitudes. Gain instabilities are corrected monitoring the time 
behaviour of thermal pulses of fixed energy, generated every 
few minutes across a Si heater resistor attached on the crystal 
absorber [23, 24]. Finally the amplitude to energy conversion 
(calibration) is determined by measuring the pulse amplitudes 
corresponding to fixed calibration lines. In the measurements 
here reported the signal of the scintillating bolometer (we will 
refer to this signal as to the heat or thermal signal) has been cal- 
ibrated on the basis of the full energy peaks visible in the spec- 
trum collected when the detector was exposed to an (external to 
the cryostat) 232 Th source. These peaks have a nominal energy 
of: 5 1 1 , 583, 9 1 1 , 968 and 26 15 keV. Below 5 1 1 keV and above 
2615 keV the energy calibration is extrapolated. However, the 
heat response for a particles is slightly different from the /3/y 
response in scintillating bolometers iTPfl . For the molybdates 
that are here reported this heat quenching factor is lower than 
few percent while for ZnSe it is a little bit higher, about 10 per- 
cent for ~6 MeV alpha particles [16|. This miscalibration of 



the a band however does not imply an appreciable change in 
the discrimination confidence level described below. 

In the case of the light signal the energy calibration is not 
needed and we therefore present its value in arbitrary units. 

Besides the amplitude, few other characteristic parameters of 
the pulse are computed by the off-line analysis. Some of them 
are: T rise and Td ecay , TVL and TVR. The rise-time (j rise ) and the 
decay-time (jdecay) are determined on the recorded raw-pulse as 
(tgo%-tio%) and (t3o%-t9o%) respectively. TVR (Test Value Right) 
and TVL (Test Value Left) are computed on the optimally fil- 
tered pulse A(t). They are the root mean square differences 
between the current signal A(t) and the reference pulse after 
OF filtering A (t)= H(t) ® S(t). In more detail, the filtered re- 
sponse function Ao(t) is synchronized with the filtered signal 
A(t), making their maxima to coincide, then the least square dif- 
ferences of the two functions are evaluated on the right (TVR) 
and left (TVL) side of the maximum on a proper time interval 
which is usually chosen depending on the shape of the OF sig- 
nals. Although these two parameters do not have a direct phys- 
ical meaning, however they are very sensitive (even in noisy 
conditions) to any difference between the shape of the analyzed 
pulse and the response function. Consequently, they are used 
either to reject fake triggered signals (e.g. spikes) or to iden- 
tify variations in the pulse shape with respect to the reference 
response function (and this will be our case). 

4. Pulse shape signature in the heat pulse 

A series of measurements was carried out, in which different 
scintillating bolometers, each coupled to a light detector (but in 
one case), were exposed to y and a sources. This allowed us 
to study the response of our devices to different radiations. The 
light signal was used to identify - on the basis of the heat to 
light ratio - the particle producing the event under study. 



As mentioned in section 3.2 for each triggered signal differ- 



ent pulse shape parameters are computed by the off-line anal- 
ysis, generally to isolate spurious and pile-up events. In the 
case of scintillating bolometers, looking at the distribution of 
the pulse shape parameters for the heat signals, we realized that 
it was possible to distinguish /3/y from a events. This is clearly 
evident in Fig. [T] where the scatter plot of the amplitudes mea- 
sured for the light and heat signals (Light vs. Heat) is compared 
with the scatter plot (obtained for the same events) of the lin- 
earized rise-time (see later in the text) vs. amplitude for the heat 
signal (t' ! " vs. Heat). In this detector, a and /3/y interactions 
draw different distributions in both the scatter plots, definitely 
proving that the shape of the thermal pulse induced by an a 
particle is different from that of a /3/y interaction. 

This behavior can be explained by the dependence of light 
yield on the nature of the interacting particle. The high ioniza- 
tion density of a particles implies that all the scintillation states 
along their path are occupied. This saturation effect does not 
occur or at least is much less for /3/y particles. Therefore, in 
a interactions a larger fraction of energy flow in the heat chan- 
nel with respect to /3/y events. This leads not only to a differ- 
ent light and heat yield but also to a different time evolution of 
both signals. The pulse shape of the thermal signal then can be 
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Figure I: CaMoC>4 crystal faced to two a sources. Left panel: scatter plot of Light vs. Heat. The energy calibration is performed on the y peaks. In red (light in the 
bottom curve) events in the 2-4 MeV region due to a particles while in green (light in the top distribution) events due to fijy in the same region. Right panel: Scatter 
plot of r? vs. Heat in CaMoC>4 crystal obtained for the same events. In green (light in the upper curve) events in the 2-4 MeV region due to fijy particles and in 
red a events. 



explained by the partition of energy in the two channels with 
different decay constants. In particular, as shown by 11281 l29l 
the scintillation produced in molybdates by a and fijy particles 
presents few decay-time constants with different relative inten- 
sities. Also a strong temperature dependence of the averaged 
decay-time of the light pulses in these same crystals was re- 
ported ll30l[3Tl . For this reason - while at room temperature the 
averaged decay-time of molybdate or tungstate scintillators is 
of the order of tens of ji/s (thus almost instantaneous in the heat 
pulse timescale) - at low temperatures it increases to hundreds 
of fi& and is therefore comparable to the typical rise-time of the 
heat signal of our scintillating bolometers. 

In the following sections we analyse the results obtained 
with different crystals, each being a possible candidate for a 
f3fJ(0v) experiment. We quote for each crystal the discrimina- 
tion power achieved - in the f3fJ(Qv) region - between a and fjjy 
particles on the basis of the light/heat ratio or simply on the ba- 
sis of the pulse shape of the heat signal. We discuss in more 
detail the case of CaMoC>4, briefly summarizing the results ob- 
tained for other crystals. 

5. CaMo0 4 

Recently CaMoC>4 has been intensively studied, for its possi- 
ble application as a scintillating bolometer for f3f3(0v) and Dark 
Matter experiments JSJ [34] [35] . This crystal contains two iso- 
topes that could undergo f3f3(0v) : 48 Ca (0^=4.27 MeV) and 
ll)0 Mo (0^=3.03 MeV). Actually, while the large content of 
ll)0 Mo makes this crystal very attractive, the presence of 48 Ca 
is a problem. Indeed, the natural isotopic abundance of 48 Ca 
(a.i.=0.19%) is too low to study the f3/3(0v) without enrichment, 
which is extremely difficult, and at the same time, it is too 
high to study the /3fJ(Qv) of I00 Mo, since the background due 
to the fJf3(2v) of 48 Ca in the f3f3(0v) region of 100 Mo will limit 




Light [a.u.] Linearized Rise Time [ms] 



Figure 2: Evaluation of the fi/y and a discrimination power in CaMoC>4. Left 
panel: projection on the Light axis of the 208 T1 y-line (2615 keV) (in green on 
the right) and of a events (in red on the left) that released a similar energy in 
the CaMoC>4 crystal. The separation corresponds to 12.6 sigma. Right panel: 
projection on the tjSL axis of events in the 2-4 MeV regions due to ply events 
(in green on the right) and a events (in red on the left). The separation here 
corresponds to 6.5 sigma. These same events are reported, with identical colors, 
in Fig.[T] 

the reachable sensitivity for the latter isotope. A possible so- 
lution to this problem was suggested by Annenkov et al. |36| 
who proposed an experiment with CaMoC>4 depleted in 48 Ca. 
Despite these possible problems, we discovered that CaMoC>4 
is an extremely interesting crystal because of its capability to 
discriminate fijy from a, thanks to the different shape of the 
thermal pulses. 

The sample we used was a cylindric CaMoC>4 crystal with 
a mass of 158 g (h = 40mm, = 35mm). The crystal was 
faced to two a sources. Source A was obtained by implantation 
of 224 Ra in an Al reflecting stripe. The shallow implantation 
depth allows to reduce to a minimum the energy released by 
the a's in the Al substrate so that monochromatic a lines (those 
produced in the decay chain of 224 Ra to the stable 208 Pb isotope) 
can be observed in the scintillating bolometer. Besides these a 
particles - all with energies above 5 MeV - the source emits a 
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/3 with a maximum energy of 5 MeV, due to the decay of 208 T1. 
Since our main goal was to study the efficiency of a particle 
rejection in the /3/3(0v) region (i.e. at about 3 MeV), a second 
source (B) was also used. This was obtained contaminating an 
Al stripe with a 238 U liquid solution and later covering the stripe 
with an alluminated Mylar film (6 fim thick). Thus the source 
produced a continuous spectrum of a particles, extending from 
about 3 MeV down to 0. 

In Fig. [T] we show the Light vs. Heat scatter plot, collected 
while the crystal was exposed to an external 232 Th source. The 
total live time of this measurement was ~43 h. The FWHM en- 
ergy resolution, on the heat signal, ranges from 2.7 keV at 243 
keV to 8.7 keV at 2615 keV on /3/y events and is about 10 keV 
on 5 MeV ar's. The light yield of the crystal, evaluated cali- 
brating the light detector with a 55 Fe source, is 1.87 keV/MeV 
(for more details on the procedure used for the evaluation of the 
light yield see references lfl4l[T6l ). 

The two separate bands - clearly visible in the Light vs. Heat 
scatter plot - are ascribed to /3/y's (upper band) and cr's (lower 
band). The upper band is dominated by the 232 Th source y's 
(plus the environmental y's). The lower band is due to a's from 
source A and from Uranium and Thorium internal contamina- 
tion of the crystal (these are the monochromatic lines above 
4 MeV) and source B (the continuum counting rate below 4 
MeV). Above 8 MeV we observe a group of events ascribed to 
a+/3 summed signals due to internal contamination in Bismuth 
and Polonium. Indeed, due to the long rise-time of this device 
(~5 ms), the beta decay of 214 Bi or 212 Bi (respectively of 238 U 
and 232 Th chains) followed immediately by a decay of 214 Po 
(t=164 fj.s) and 212 Po (r=298 ^s ), may lead to a pile up on the 
rise-time of the thermal pulses that can hardly be recognized 
as a double signal. The two decays produce therefore a single 
pulse, with an energy that is the sum of the two. 

The discrimination between the a and the /3/y populations, 
provided by the scintillation signal, can be evaluated by mea- 
suring the difference between the average amplitude of the light 
signal produced by the two kinds of particles (Lightg/ y and 
Light„), considering a group of events releasing the same en- 
ergy in the scintillating crystal. This difference is then com- 
pared with the width of the two distributions (o~p/ y and cr a ). 
The discrimination confidence level Du g ht can be then defined 
as: 



D 



Light 



Lightpiy - Light a 



4 



(4) 



To evaluate this discrimination power in the /3/3(0v) region we 
selected events belonging to the 2615 keV y-line full energy 
peak and compared their light pulse distribution with that of 
events of similar energy in the a band (Fig. |2]left panel). T>ught 
results to be 12.6 sigma. The reason for using events with the 
same energy to evaluate D/,^, is in the energy dependence of 
the distance (and width) of the a and /3/y bands, that induce an 
energy dependence of the discrimination confidence level. In 
the energy range where both the bands are populated (i.e. be- 
tween 1 and 3 MeV), D L , g /„(E) appears to be linearly decreasing 
with energy. Its extrapolated value at E=0 being 3 sigma. 




Figure 3: CaMoCXi crystal, comparison of /3/y's and a's pulse shape. Top panel: 
thermal pulse obtained by averaging signals belonging to the 2615 keV y-line 
(black continuous line) and signals due to a particles that release a similar en- 
ergy (red dotted line). Bottom panel: difference, in percent of the maximum 
amplitude, between the two averaged pulses. 



As already anticipated, we discovered that the heat signal 
shape is enough to discriminate /3/y's from a's. This is shown 
in the right panel of Fig. [T]where we report the (thermal pulse) 
t' ! " vs. Heat scatter plot for the same events shown in the left 

rise 1 

panel. Two separate bands, ascribed to /3/y and a events can be 
identified. The former with an average rise-time of ~5.8 ms, 
the latter with ~5.6 ms. In this plot the weak energy (Heat) 
dependence of r rise observed for the two populations was cor- 
rected by fitting their rise-time distributions with two lines hav- 
ing the same slope. To do this we used /3/y events in the 0.5-2.6 
MeV range and a events in the 1.5-7 MeV range. The T nje is 
then linearized re-defining it as t5J = t>,. sc - slope x E with 
slope=0.0075 ms/MeV. 

In order to emphasize the correspondence in the identifica- 
tion of /3/y and a events, signals selected in the Light vs. Heat 
scatter plot as having an energy between 2 and 4 MeV and be- 
longing to the /3/y or a bands, are reported in different colors. 
We will refer in the following to these two groups of events as to 
2-4 MeV y's and 2-4 MeV a's, although - above the 2615 keV 
line - the y group is empty. These two samples will be used to 
evaluate the difference in shape among /3/y and a signals. 

In Fig. [3] we compare two pulses obtained by averaging - sep- 
arately - signals belonging to the 2615 keV y-line and signals 
due to a particles that release a similar energy. The average 
is here needed to get rid of the noise that can mask the small 
differences among the pulses. The difference in shape can be 
appreciated in the bottom panel of Fig. [3] 

The discriminating power T>Ri se Time of the heat pulse shape 
method can be defined exactly as done for T>u g ht'- 



D 



RiseTime 



(E) 



RiseTimepiy — RiseTime a 



4 



(5) 



where, however, the dependence on energy is here due only to 
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Figure 4: Z11M0O4 crystal. Scatter plot of Light vs. Heat (left panel) and scatter plot of RDC vs. Heat (right panel). In green events in the 2-4 MeV regions due to 
Ply particles and in red those due to ci's. 



the variation of the width of the rise-time distributions since 
the difference in rise-time appears to be independent on energy. 
The distribution obtained projecting the rise-time of 2-4 MeV 
y's and 2-4 MeV a's is shown in Fig. [2] (right panel). A gaus- 
sian fit of the two peaks yields a rise-time of (5.788±0.017) ms 
for y's and of (5.649+0.013) ms for a's. T>Ris e Time results equal 
to 6.5 sigma. As for the T>ugkt case, we can extrapolate the 
D RiseTime value also in the energy region where the a band is 
poorly populated, simply assuming that the distance between 
the two bands remains constant (as we observe) and Dja se Ti me 
changes because the width of the two distributions becomes 
larger, at low energies, due to noise. The result is that Y^RiseTime 
becomes lower than 2 sigma at 500 keV. In other words the 
~150 fis difference in rise-time of a and fi/y events cannot be 
appreciated when the pulse amplitude is too low and, therefore, 
the noise modifies appreciably the signal shape. Finally, al- 
though in this measurement the rise-time is the most efficient 
shape parameter for a event discrimination, good discrimina- 
tion levels were observed also in other parameters such as the 
decay-time and the TVR. The possibility to increase the dis- 
crimination power by combining the information from different 
shape parameters or the fit on the rise-time is under study. 

We note that either in the Light vs. Heat scatter plot as in 
the T^" se vs. Heat one we can see some outliers that could in- 
dicate a possible failure of the particle identification technique. 
However such events can be accounted for if one considers the 
following two effects. If a y-ray interacts both in the light de- 
tector and in the scintillator the light signal that is read-out has 
a wrong amplitude since it is not only ascribed to scintillating 
photons but also to a direct y interaction in the Ge wafer. The 
measurement is performed in a high rate condition, therefore 
we expect to have a number of pile-up events in each of the 
two detectors, this leads to an erroneous evaluation of the pulse 
amplitude and pulse rise-time. 

To conclude, in the case of /3y8(0v) decay the discrimination 
power provided by the use of the scintillation signal is compa- 



rable with that provided by the pulse shape analysis. For what 
concerns the use of this device for the measurements of a par- 
ticle emissions from an external sample, a D RiseTime better than 
2 sigma above 500 keV means that this detector has a good 
sensitivity even in the region where Si surface barrier detectors 
start to be dominated by y background. On the contrary the ap- 
plicability of this technique to Dark Matter searches cannot be 
proved directly. 

5.1. Other crystals 

Other scintillating crystals have shown the possibility to 
discriminate interacting particles through the thermal pulse 
shape differences. Among the tested crystals, other molyb- 
dates (Z11M0O4, MgMo04) and also other crystals such as ZnSe 
showed a good discrimination power. 

5.2. ZnMoO A 

We tested a 19.8 g ZnMo04 crystal, having the shape of a 
prism with height of 11mm and a regular hexagonal base, with 
a diagonal of 25 mm 11171 . The FWHM energy resolution, on 
the heat signal, is 4.2 keV FWHM on Ply events at 2615 keV, 
and about 6 keV on the 5.4 MeV a line. The light yield of the 
crystal is 1.1 keV/MeV. The total live time of this measurement 
was ~ 195 h. Also for this analysis we have used the Light vs. 
Heat scatter plot to select fi/y and a between 2 and 4 MeV, in 
order to evaluate the discrimination power provided by pulse 
shape studies. 

Unlike for the CaMo04 case where the most efficient param- 
eter for the discrimination is the rise-time, this measurement 
showed a higher discrimination power on the decay-time of the 
thermal pulse. In order to emphasize the discrimination power, 
all signals have been fitted with a function obtained as the sum 
of two exponentials: 

AV(f) = ( e - ,/Ti+A < + e - t/Ti+M ) (6) 
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Figure 5: Evaluation of the fj/y and a discrimination power in ZnMoQ^ Pro- 
jection on the RDC'" 1 axis (after linearization of the scatter plot of Fig. Ml right 
panel) of events in the 2-4 MeV region due to ft/y events (in green) and a events 
(in red). 




Linearized Ratio Decay Constants 



Figure 6: MgMo04 crystal: scatter plot of RDC vs. Heat (top panel) and 
projection of 2-4 MeV events on the RDC'"' axis (bottom panel). In green 
events in the 2-4 MeV region due to fijy particles and in red those due to a's. 



where the ti and T2 parameters are obtained by fitting the decay 
of the thermal signals in the raw-pulse acquired for each event. 
It was observed that the best discrimination power is obtained 
by the ratio of the two decay constants (RDC): 

RDC = — (7) 

The scatter plot of RDC vs. Heat is reported in the right panel of 
Fig. |4](in the left panel the corresponding Light vs. Heat scat- 
ter plot). The discrimination between a and fijy populations 
provided by RDC is evaluated linearizing the RDC vs. Heat 
relationship (fitting the a band in scatter plot of Fig. |4]with a 
polynomial as done for the r nje of CaMoCu) and then project- 
ing the distribution of RDC' m for events in the 2-4 MeV range. 
The two peaks (fijy and a) are then fitted with a gaussian (Fig. 
|5j. The discrimination power Drdc is 6.4 sigma. 

5.3. MgMo0 4 

This compound contains, as in the two previous cases, the 
f3f3 active isotope l00 Mo which is here present in a larger con- 
centration (~52% in mass). 

The crystal tested is 32x31x24 mm 3 with a weight of 89.1 
g. The total live time of this measurements was ~22 h. In this 
run it was not possible to face the crystal to a light detector be- 
cause of the assembly structure so we can't use the Light vs. 
Heat scatter plot in order to tag a events. The performances of 
the bolometer were quite poor, most probably this was due to a 
problem with the gluing of the NTD thermistor: at the end of 
the measurement, when the crystal was back to room tempera- 
ture, we discovered that under the thermistor the crystal showed 
a crack. This could explain why the signal to noise ratio was so 
bad (the energy resolution measured on a lines was ~150 keV 
FWHM) and consequently also the resolution in the evaluation 
of the signal shape parameters was limited. Despite these prob- 
lems, we were able to observe in the scatter plot of RDC'"' vs. 



Heat (fig. [6]) clear difference between events due to the interac- 
tion of a particles and events due to fijy events. 

The discrimination between the two populations provided by 
the RDC parameter can be evaluated by means of a gaussian 
fit of the two peaks obtained by projecting selected events after 
linearization. The resulting discrimination power Drdc is 1-8 
sigma. 

These preliminary results, even if limited because of the 
problems reported above, lead us to program new more detailed 
measurements to better study this promising crystal. 

5.4. ZnSe 

&2 Se is a fif! emitter with an isotopic abundance of 9.2% and 
a Q-value of (2995.5 ± 2.7) keV. It has always been considered 
a good candidate for f3f3(Qv) studies because of its high transi- 
tion energy and the favorable nuclear factor of merit. For these 
reasons in the last years an R&D work was carried out in which 
we have extensively studied the performances of ZnSe detectors 
in different conditions. For our studies we have used different 
ZnSe crystals. Characteristics of measurements done and ob- 
tained results are reported in detail in [ 1 6j . Here we report the 
discrimination power on the RDC obtained with the Huge ZnSe 
crystal (h = 50mm, = 40mm, 337g). 

In Fig. [7]the Light vs. Heat scatter plot of a measurement of 
~70 h of live time is shown. In order to have a high number 
of a counts in the 2-3 MeV region also in these measurements 
a degraded 238 U source was placed in front of the crystal. An 
explanation of the Quenching Factor larger than one (i.e. the a 
band lies above the fijy band) and other anomalies observed in 
this crystal can be found in |[T6l . 

Also in this case, the RDC'"' vs. Heat (fig. [8} showed a 
difference between events due to the interaction of a particles 
and events due to fijy events. The discrimination power Droc 
is 2.2 sigma. 
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Figure 7: ZnSe Crystal: scatter plot of Light vs. Heat. In green events in the 
2-4 MeV region due to ft/y particles and in red those due to a's. 

6. Conclusion 

The possibility to discriminate the nature of the particle inter- 
acting in a bolometric detector, simply on the base of the shape 
of the thermal pulse, is now definitely proved and opens new 
possibilities for the application of these devices in the field of 
rare events searches. In particular, the high rejection capabil- 
ity that could allow to completely rule out the a background in 
j6/S(0v) experiments was demonstrated. Unfortunately, based on 
the results obtained so far, the applicability of this technique to 
Dark Matter searches cannot yet be proved directly. This fea- 
ture was observed in different scintillating crystals (CaMoC>4, 
Z11M0O4, MgMo04 and ZnSe) and new tests are under prepara- 
tion in order to investigate if a similar behavior can be observed 
also in other compounds. A discrimination confidence level 
of ~6.5 sigma was obtained both for CaMoCU and ZnMoC>4 
crystals in the 2-4 MeV energy region. Discrimination con- 
fidence levels reached with MgMoC>4 (D SDC =1.8cr) and ZnSe 
(D RDC -2.2(f) are not so high but they indicate that even with 
these crystals the discrimination based on pulse shape analysis 
is possible. New techniques aiming at improving the discrimi- 
nation power of the pulse shape analysis are being studied. 
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